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We show in a mean-field calculation that phase diagrams remarkably similar to those recently 
proposed for the cuprates arise in simple microscopic models of interacting electrons near half- 
filling. The models are extended Hubbard models with nearest neighbor interaction and correlated 
hopping. The underdoped region of the phase diagram features d^2_y2 density-wave (DDW) order. 
In a certain regime of temperature and doping, DDW order coexists with antiferromagnetic (AF) 
order. For larger doping, it coexists with d^i_y2 superconductivity (DSC). While phase diagrams of 
this form are robust, they are not inevitable. For other reasonable values of the coupling constants, 
drastically different phase diagrams are obtained. We comment on implications for the cuprates. 

PACS numbers; 74.72.-h, 71.10.Hf, 74.25.Dw 



I. INTRODUCTION 

The high- Tc cuprates exhibit peculiar behavior when 
underdoped: the density-of-states is depleted at low en- 
ergies, as if some of the degrees of freedom of the system 
were developing -a gap. STbis behavior, observed in optk 
cal conductivity ,El NMR,H angle-resolved photoemission,Ll 
c-axis tunnelingjQ and specific heat measurements Jj was 
dubbed the 'pseudogap'. The emergence of the 'pseu- 
dogap' mimics somewhat the impoverishment of the low- 
energy excitation spectrum which accompanies the devel- 
opment of dr^2_y2 superconductivity and resembles, more 
generally, the type of gap formation which is concomi- 
tant with a large class of order parameters. However, 
it does not - at first glance - appear to be connected 
with the formation of an ordered state. Consequently, it 
was initially believed that the pseudogap was a crossover 
phenomenon, and attempts to describe it depended on 
various approximate methods of treating states with lo- 
cal, fluctuating order.Ql3 

However, it has recently been proposed that the 'pseu- 
dogap' state is actually a broken-symmetry ordered 
state, and that the signatures of the order are s 
enough that the state was able to appear incogniti 




In Ref. H, the dx2-y2 density- wave (DDW) stateEj was 
advanced as a candidate order. The realization that this 
is a realistic possibility has led to a re-examination of 
the experimental circumstances. Recent elastic neutron 
scattering experiments, which directly probe the symme- 
tries broken by DDW order - time-reversal and trans- 
lation by one lattice spacing - appear to have observed 
it 113 A number of other experiments are consistent with 
the proposal,El especially measurements of the superfluid 
density as a function of doping £3 

The experimental situation seems promising, which is 
strong incentive to reconsider the theoretical state of af- 
fairs. If the 'pseudogap' state is, indeed, an ordered state, 
then we should be able to study it within mean-field the- 
ory, as we would study the antiferromagnetic state, su- 



perconducting state, or any other ordered state. Mean- 
field theory is unlikely to explain the detailed shape of the 
phase boundary, but one can hope that it will capture the 
broad features of the phase diagram, such as its topol- 
ogy and the basic temperature scales. Deep within any 
phase, with T ^ and far from any quantum phase tran- 
sitions, the mean-field Hamiltonian should be the correct 
Hamiltonian, although the parameters in it may need to 
be renormalized from their mean-field values. Thus it 
seems natural to simultaneously study the antiferromag- 
netic (AF), dx2-y2 -wave superconducting (DSC), and 
DDW order parameters in mean-field theory. The inter- 
play and possible coexistence of these orders should be 
qualitatively and semi-quantitatively explained by mean- 
field theory. Phase transitions, quantum or thermal, may 
not be accurately described in their asymptotic limits, 
but the AF, DDW, and DSC phases will, as will possible 
phases with coexisting AF, DDW, and DSC orders. 

However, there is an immediate problem faced by such 
a program. What microscopic Hamiltonian should be 
used? In the early days of high- Tc, it was hoped that the 
important physics of strong local repulsion and superex- 
change, which is present in the simplest models, such as 
the Hubbard and t — J models, would be sufficient to ex- 
plain all of the interesting physics of the cuprates. This 
appears not to be the case. Monte Carlo studies haiie 
not found superconductivity in the Hubbard model,Ej 
while Monte Carlo calculations, exact diagonalization, 
and density-matrix renormalization-group (DMRG) cak 
culations give conflicting results for the t — J model.ta 
DMRG studies have found that the behavior of n-leg 
ladders depends sensitively on iJie strength of, for in- 
stance, |Sacond-neighbor hopping,E3 as have Monte Carlo 
studies.EZI Indeed, some numerical poults are sensitively 
dependent on boundary conditionsJla which is further in- 
dication of the instability of many of these models to rel- 
atively small changes in the parameters. Furthermore, 
the physics of charge-ordering is probably not correctly 
described by the t — J model without near- neighbor (and 
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possibly long-range) Coulomb repulsion.tSESO Indeed, 
it is also clear from experiments that relatively small 
changeg-p such as those associated with substituting Nd 
for La,E3 which is off the radar screen of the t — J and 
Hubbard models - can radically change at least some as- 
pects of the behavior of these materials. In short, the 
detailed form of the underlying Hamiltonian matters. 



Fortunately, we are not completely in the dark about 
the nature of the microscopic Hamiltonian. Local 
Coulomb repulsion, both on-site and near-neighbor, is 
clearly an important part of the physics. This is known 
from microscopic calculations of the Hubbard parame- 
ters t, U, and also from the fact that the undoped parent 
compounds are antiferromagnetic insulators. The other 
important clue - which derives entirely from experiments 
- is that the cuprates superconduct. The correct micro- 
scopic model (or models) must support d-wave supercon- 
ductivity when doped away from half-filling. If the Hub- 
bard and t — J do not - and it appears that they do not 
for t/U small - then they cannot describe the cuprates 
fully. 



Our strategy will be to take a generalization of the 
Hubbard model which includes next-neighbor repulsion 
and, most importantly, pair-hopping (or correlated hop- 
ping). The pair-hopping term favors superconductivity. 
Even when it is relatively small, it stabilizes supercon- 
ductivity in the Hubbard model, as we will see. There 
are a variety of ways in which such a term - or another 
term with similap-effect - could arise, either from quan- 
tum chemistryME^ or in the passage to an effective de- 
scription such as the t — J model; in both cases, it is 
essentially a result of strong local Coulomb repulsion, 
as superexchange is. In any event, it appears that such 
physics is necessary to stabilize superconductivity, so we 
will incorporate it in our model. We will find that such 
a term also leads to DDW order. 



To summarize, we consider a model which is chosen 
so that it incorporates the basic physics of strong local 
repulsion and so that will have a phase diagram which 
includes AF at half-filling and DSC at some finite doping. 
We find that it naturally supports DDW order. In mean- 
field theory, we find a phase diagram in the temperature- 
doping plane which resembles the experimental phase di- 
agram of the cuprates, with the DDW phase boundary 
playing the role of the experimental pscudogap onset line. 
This DDW line continues into the DSC state, so that 
the underdoped superconducting state is characterized 
by both DSC and DDW orders. At low doping, there is 
also a region of coexistence between AF and DDW or- 
ders. We comment on the interpretation of experiments 
vis-a-vis our findings. 



II. MODEL HAMILTONIAN 

We consider the pbllowing bilayer lattice model of in- 
teracting electronsO 



H — Hkin+'Hi, 



(1) 



where 

"Hkin 



(A) 



h.C, 
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^(2) , Jl)t (2) 
a ia ' i~{-x—y,a i(T 



„(l)t (2) (l)t (2) 



h2S,CT 



4a + x^y+1^2 + h.C. 



(2) 



and 



JA)t JA) (A)t JA) 



^icr" ^ja ^ja ^i'a' (3) 



In the formulas above is hopping with = t for 
nearest neighbors, = t' for next nearest neighbors and 
tij = otherwise. The other parameters are the tun- 
neling t^, the on-site repulsion U , the nearest-neighbor 
repulsion V ^ and next-nearest-neighbor correlated hop- 
ping tc- The indices i, j correspond to a lattice site, a to 
the spin, and A to the layer. 

The next-nearest-neighbor correlated hopping term is 
physically kinetic, but since it is also quartic, we are go- 
ing to treat it as interaction. It hops an electron from 
i' to j when j is vacated by an electron hopping to i. 
These two hops are correlated by virtue of Coulomb in- 
teraction between the electrons. The presence of this 
term in the pjiprates has been shown in band-structure 
calculations .l3 Correlated hopping has been discussed in 
Refs. |4|j26|j2^,|| as a possible mechanism of supercon- 
ductivity, but it has also been foundnJ that it favors DDW 
order as well. v-\t-\ 

The tunneling term is momentum conserving. E3'l3 We 
consider a Cu02 bilayer because the pseudogap has been 
best characterized in bilayer materials such as YBCO and 
Bi2212. 

To derive a mean-field theory, it is convenient to take 
the Fourier transform of Eq. (|l]) and regroup the terms. 
This task would be particularly simple if there were only 
one phase at a given set of parameters. For example, a 
DDW reduced Hamiltonian would look like 



(4) 



fc,fc' 



where / (fc) = cos fc^, — cos ky (the lattice spacing has been 
set to unit) and the DDW mean-field coupling constant 
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5DDW = 8V + 2Atc 



(5a) 



Similar values of the mean filed coupling constants can 
be derived for other phases as well. Thus, for antifer- 
romagnesm, d-wave superconductivity and (tt, tt) charge- 
density wave we derive: 



5AF 
ffDSC 
5CDW 



2U, 



2U. 



(5b) 
(5c) 
(5d) 



In fact, the interaction part of the Hamiltonian Eq. (|^) 
can be further generalized to include the interlayer 
Coulomb interactions: 



V 



'Y.-^'nf\ (6) 



where X ^ X'. Then for the given interlayer configura- 
tion of the order parameters (antisymmetric for AF and 
DDW and symmetric for DSC), the mean field coupling 
constants become: 



5DDW = 8V + 8V' + 24i„ 

ffAF = 2(7 + 2(7', 
gusc = 12U~&V + 8V'. 



(7a) 
(7b) 
(7c) 



For the opposite configuration (symmetric for AF and 
DDW and antisymmetric for DSC) the contributions of 
U' , V' would be negative, which is the main reason why 
such configurations have generally higher energy and are 
not observed. On the other hand, the fact that five in- 
teraction terms produce only three phases means that we 
can have the same phase diagrams (corresponding to a 
given set of ffp's) for a range of values of the interaction 
constants. In the following section we will assume that 
U' ^ V = so that each phase diagram will correspond 
to a unique set of U, V, tc- 

The total Hamiltonian contains the reduced parts cor- 
responding to these phases as well as the interactions 
between the order parameters. However, since we expect 
5CDW to be negative so that the corresponding order pa- 
rameter is always zero, we will ignore the term corre- 
sponding to this phase. The final form of the reduced 
Hamiltonian is 



cd = ' 
Jk 



n-i I (^)t ) 

Ik 



i/AF ; ^k+Q,cr^ka ^k'cr'^k'+Q,. 
J k,k' 

f /(fc)/(fc')4^L4t^4'^'4'i 



k,k' 



VQ.cy' 



-gusc f{k)nk')c'^^^'c'^^l\Ay%,, (8) 

J k,k' 

where Ckii = efe22 = Cfe = -2t (cos fc^; -I- cos fc^^) - 
4t' cos cos ky, efei2 = efc2i = efcj. = (t±/4)/(fc)^. 



The standard Hubbard-Stratonovich mean-field- 
theoretical treatment of Eq. (||) is to assume the 
presence of a bosonic mean field, defined as an order 
parameter, neglect the fiuctuations, find the eigenvalues 
of the Hamiltonian and finally, integrate out the fermion 
degrees of freedom to derive the free energy. 

We define the order parameters of DDW, AF and DSC 
phases as follows: 



J k 



Aa = .9DSC / / (fc) c 



„(A)tJA)t 
-kl- 



(9a) 
(9b) 
(9c) 



We assume that and M\ are anti-symmetric in the 
bilayer index. Then, the free energy of the system is 



lAfl 



lAI 



ffAF 5DDW ffDSC 



/fex>0 

si,S2,S3 = ±l ky>k^ 



2rin(^2cosh|^ [{/(fc)A}'+ (si{[efe + S2efcJ 



9-1 1/2 \' 
+ [fik)cp + s;Mfj -mJ 



1/2- 



(10) 



As we expand this expression for small values of the 
order parameters, we can construct a Landau-Ginzburg 
theory: 

/ (T) = /o (T) + ^ ap \<i>/ + J2bp \<^pf 
p p 

where p denotes the kind of the order parameter (AF, 
DDW or DSC) and is the order parameter (M, cj) or 
A, respectively). The Up coefficients cross zero at the 
transitions so that Qp — are the equations that deter- 
mine critical temperature T^. 



Sl,S2,S3 = ±l ky>k^ 



Kp (k) , (12) 



where 



i^DDW (fc) = /(fc)'i^AF (fc), 



(13a) 
(13b) 



i^DSc(fc) = :,^?^tanhf^i%^). (13c) 



2 \siek - fJ.\ 



2T 
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Here = + S2£k±- The bp coefficients are positive: 



bp = 



E 



fcx>0 



where 



Kowik)^f{kfK'^^ (fc) 



^DSC (k) — 



fikfm (k) 



\siek - Ml 



|3 ' 



where 



rji (fc) ~ tanh 
772 (fc) = tanh 



2T 



\^k\/2T 



cosh (^i^) 
Isiefc - Ml /2r 



cosh 



.h(^i^)' 



(14) 

(15a) 
(15b) 
(15c) 

(16a) 
(16b) 



Finally, the c„„' coefficients are 



Cppi — 

Sl,S2,S3 = ±l ky>k. 



k^>0 ^PP' 



where 



-f'^AF.DSC (k) 



fiky 



A\slek-^J.\ \ek 

^''^DDW.DSC ^ / (^) ^AF.DSC (^) 
^AF.DDW (k) = ^]Z-^m {k) . 

4 kfc 



?72 (fc) , 



(17) 

(18a) 
(18b) 
(18c) 



The fact that all K^p, > implies that the phases com- 
pete with each other. 



III. PHASE DIAGRAM 

The mean field phase diagram can be derived by min- 
imizing Eq. ( |lO| ) at fixed doping. 

Since there is a large number of parameters in our 
model, there is substantial variety in the possible dia- 
grams. One such diagram, generated with t = 0.5 eV, 
t' = -0.025 eV, ti_ = 0.05 eV, 5af = 0.06 eV, c^ddw = 
0.02 eV, 5DSC — 0.01 eV, is shown on the figure ||. 
The corresponding values of the interaction constants are 
U ~ 0.03 eV, y = eV, ~ 0.8 x lO^^ eV. Note that 
for these values of the constants, 5cdw = —0.04 eV < 0, 
which is consistent with our assumption that a (tt, tt) 
charge-density wave is not energetically favorable. 

Another diagram, shown on the figure H, was generated 
with t = 0.5 eV, t' = 0, t± = 0.1 eV, ^af = 0.084 eV, 
5DDW = 0.038 eV, and ^dsc = 0.017 eV. The correspond- 
ing values of the interaction constants are U ~ 0.042 




0.1 0.15 0.2 
doping 

FIG. 1: Phase diagram for t = 0.5 eV, t' = -0.025, tj_ 
0.05 eV, U ~ 0.03 eV, F = eV, tc ~ 0.8 x 10"^ eV. (^af 
0.06 eV, goDw = 0.02 eV, ^dsc = 0.01 eV.) 




0.15 0.2 

doping 

FIG. 2: Phase diagram for t = 0.5 eV, t' = 0, fx = 0.1 eV, 
U ~ 0.042 eV, V ~ 1.7 X 10"'' eV, ~ 1.5 x 10-=* eV. 
(3AF = 0.084 eV, c^ddw = 0.038 eV, ^dsc = 0.017 eV.) 



eV, ~ 1.7 X 10-4 eV, ~ 1.5 x lO'^ eV, and also 
gcBW - -0.045 eV. 

As we can see, in both diagrams the antiferromag- 
netic transition temperature at half-ffiling is close to 1000 
K. This should be understood as the scale at which 
two-dimensional antiferromagnetic correlations develop 
locally. Due to the Mermin-Wagner-Coleman theorem, 
which states that a continuous symmetry cannot be bro- 
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ken spontaneously at finite-temperature in 2D, the tran- 
sition temperature is zero for a single bilayer. The cou- 
pling between different bilayers (which is not included 
in our single-bilayer calculation) stabilizes the antiferro- 
magnetic phase with a transition temperature around 410 
K. In lightly-doped cuprates, the presence of impurities 
causes the misalignment of locally ordered antiferromag- 
netic clusters, thereby forming a spin glass. Thus, if we 
interpret our T/v as the scale of local 2D antiferromag- 
netic order, which could become 3D antiferromagnetic 
order or spin glass order, then the phase diagrams of 
Figs. |l| and ^ are very reasonable, indeed. 

Experiments might lead us to expect that DDW order 
would occur in the range of doping between 0.07 and 
0.19. This range is smaller than one shown on Fig. |^ and 
a bit larger than that shown in Fig. |l|. The temperature 
scale for this phase on Fig. |l| is very reasonable; it is 
almost three times higher on Fig. ||. This change occurred 
primarily as a result of the increased value of tc- If we 
further increase tc to 1.9 x lO^'^ eV, the DDW phase 
will begin to suppress the AF phase and will expand up 
to half-filling at finite temperatures. In general, varying 
the interaction constants by less that 20 - 30% does not 
change the phase diagram qualitatively. However, larger 
variations lead to completely different classes of phase 
diagrams, such as those with the AF phase suppressed or 
without a DDW phase at all. For example. Fig. || shows 
the case when due to the smaller value of correlation 
hopping, both DDW and DSC phases disappear and only 
AF phase remains in the diagram. 
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AF 
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0.15 



0.2 



0.25 



0.3 



doping 



FIG. 3: Phase diagram for t = 0.5 eV, t' = -0.025, ti_ 
0.05 eV, U = 0.05 eV, = eV, U =i 3.3 x 10"* eV. (^af 
0.1 eV, gDDw = 0.008 eV, ffDSC = 0.004 eV.) 



The DSC phase occupies a doping range away from 
half-filling primarily as a result of band structure effects 
associated with the bilayer splitting. In the absence of 
other orders, it would extend all the way to half-filling, 
but it is suppressed at low doping by DDW and AF or- 
der. In a more realistic calculation, superconductivity 
would be suppressed close to half-filling by no-double- 
occupancy constraint, i.e. by strong local Coulomb re- 
pulsion. However, the DSC phase never even makes it 
that close to half-filling because the DDW phase inter- 
venes. 

An important feature common to both diagrams is the 
existence of regions with two simultaneous kinds of or- 
der. Namely, there is a region with DDW-I-AF order and 
a region with DDW-I-DSC order. The system is an insu- 
lator in the AF state at half-filling, a metal in the DDW 
and DDW-I-AF states, and a superconductor in the DSC 
and DDW-I-DSC states. 

All of the transitions are of second order at the mean- 
field level because of the signs of the Cppi couplings be- 
tween the order parameters in the Ginzburg-Landau the- 
ory Eq. 

The calculated dependence of the chemical potential 
/i on the doping x inside the DDW phase and in its 
proximity is nonmonotonic. This is due to the rapid 
development of the DDW gap, which causes the chemi- 
cal potential to be lower than in the normal state. The 
thermodynamic inequality {dfi/dx)T,v £ implies that 
when this is violated, mean-field theory should be cor- 
rected using Maxwell's construction, which signals that 
fiuctuations drive the transitions first-order as a function 
of /Lt. Consequently, we would expect the underdoped 
side of the DSC phase to be characterized by a smjaller 
than expected chemical shift, as has been observed.L3 A 
first-order phase transition as a function of chemical po- 
tential is manifested as phase separation in a two-phase 
coexistence region spanning a range of dopings when the 
doping is held fixed instead. It has been argued that such 
phase separation will be precluded by Coixkmb interac- 
tions, thereby leading to stripe formation.EScj 



IV. CONCLUSION 

We have studied the phase diagram of a bilayer lattice 
model using mean field theory. Since we have focused on 
ordered phases, this should be a valid approximation. We 
found that for certain ranges of the values of the inter- 
action constants the phase diagram agrees well with the 
experimentally observed phase diagram of YBCO if the 
'pseudogap' is associated with DDW order. The diagram 
remains in qualitative agreement with the experimental 
data when the parameters of our model vary by less than 
20 - 30% and becomes qualitatively different for larger 
variations. Clearly, such a phase diagram is reasonably 
robust, but is hardly inevitable. This is reassuring be- 
cause high-temperature superconductivity is stable, but 
only appears in a special class of materials (to the best 
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of our current knowledge) . 

There are some systematic errors associated with 
mean- field theory, on which we now comment. It under- 
estimates the effect of fluctuations. Thus, the Neel tem- 
perature is very large in mean-field theory, while it should 
actually be zero in any strictly two-dimensional system. 
However, the Neel temperature which we find should be 
regarded as the temperature-below which a renormalized 
classical description is valid. U The Neel temperature ob- 
served in experiments is associated with the crossover 
from 2D to 3D. Mean- field theory also overestimates the 
coupling which drives antiferromagnetism, which it takes 
to be essentially U. For small U, this is correct, but 
for large U, it should be replaced by J ~ t'^/U. Indeed, 
the largc-C/ limit is generally somewhat problematic near 
half-filling since the Gutwiller constraint is not enforced 
in mean-field theory. The €1^2 _y2 symmetry of the DDW 
and DSC states lead one to the erroneous conclusion that 
they are completely unaffected by large U. This cannot, 
of course, really be true; clearly, mean-field theory under- 
estimates the tendency of large-f/ to push these ordered 
states away from half-filling. The seemingly small value 
of U taken in our calculation should be interpreted in 
light of these observations. Other mean-field treatments 
which incoporate strong local Coulomb repulsion more 
prominently have alsoJQund DDW order in a generaliza- 
tion of the t — J modelE3 andJn the Hubbard model with 
nearest neighbor attractionEj 

We find that the scale associated with superconduc- 
tivity is largely determined by the strength of correlated 
hopping. At the moment, this is rather ad hoc, but we 



had little choice but to introduce some term of this sort in 
order to have a phase diagram which includes supercon- 
ductivity. It is possible that the superexchange coupling 
J plays a more important role than we have accorded it 
in setting Tc, but superexchange is beyond a mean-field 
treatment. 

As we have seen, the very term which stabilizes su- 
perconductivity also supports the development of DDW 
order. One way of interpreting our results begins with 
the observation that the DDW order parameter, when 
combined with the real and imaginary parts of the DSC 
order parameter ripjem a triplet under an SU (2) group of 
transformations. u'Eil If this 'pseudospin' SU (2) is a sym- 
metry of the Hamiltonian, then DDW and DSC orders 
will be equally favored. Thus, one can envision that the 
important order-producing term in the Hamiltonian is 
5L/(2)-symmetric while small symmetry-breaking terms 
drive the system into either the DDW or DSC states. Our 
result shows that pair-hopping is of this form. Are all 
physically reasonable mechanisms for d^2_j^2 supercon- 
ductivity similarly invariant under pseudospin SU(2)7 
This is an open problem; we have answered in the af- 
firmative for one particular class of Hamiltonians. 
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